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ABSTRACT 


Sensitivity  of  designed  plans  and  ocrputed  risks  to  changes 
in  the  prior  parameters  has  been  studied  in  this  report.  In  the 
first  part,  the  study  is  pursued  for  specified  Changes  in  the  prior 
mean  and  the  shape  parameter.  The  second  part  involves  an  arpirical 
quantification  of  the  uncertainty  in  parameter  estimates  by  using 
Monte  Carlo  simulation,  likelihood  contours  and  the  asynptotic 
distribution  of  the  parametric  estimators. 
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1.  INTRODUCTION 


Various  aspects  of  the  design  and  analysis  of  single  sanple  reliability 
acceptance  sanpling  plans  in  the  presence  of  a prior  distribution  have  been 
described  in  Goel  and  Jbglekar  (1976  a,b,c,d) . These  reports  primarily 
deal  with  the  case  when  the  failure  distribution  is  exponential,  given  by 
f (t 1 6)  - i t » 0,  0 > 0 (1) 

where  t is  the  time  to  failure  and  0 is  the  mean  time  between  failures 
(MTBF) . The  parameter  0 is  construed  to  be  a random  variable  with  an 
inverted  gamma  prior  density 

e-^.Ky.  e>o.  <2, 

In  the  presence  of  a prior  dmnsity,  the  producer's  and  the  consoler's  risks 
can  be  quantified  in  different  ways.  The  definitions,  interpretations, 
interrelationships  and  the  appropriateness  of  various  risks  have  been 
discussed  in  Goel  and  Joglakar  (1976  b) . Determination  of  t te  prior 
density  g(0)  and  the  estimation  of  the  parameters  y and  \ have  been  des- 
cribed in  Goel  and  Joglakar  (1976  c) . Nirarical  ard  graphical  procedures  for 
the  design  of  single  sanple  plans  (T,r*)  are  discussed  in  Goel  aid  Jbglekar 
^^76  > where  T is  the  test  time  and  r*  is  the  n! Inwnhlf  number  of 

failures  in  time  T. 

ttie  purpose  of  this  report  is  to  investigate  the  sensitivity  of  the 
designed  plans  and  the  computed  risks  to  changes  in  the  prior  parameters. 

This  investigation  has  been  pursued  in  two  parts.  In  the  first  part  we 

study  the  sensitivity  for  specified  changes  in  the  prior  mean  u ( - JL_) . 

X P 


A 
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! 
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and  the  shape  parameter  A.  The  effect  of  changes  in  T*  ( = T/6Q) , r*, 
and  ypOn  various  risks  is  considered  in  Section  2.1.  The  following 
sensitivity  analyses  are  performed  in  Section  2.2: 

(i)  For  fixed  risks,  the  effect  of  small  changes  in  Ppand  A on  the 
designed  values  of  T and  r*. 

(ii)  For  plans  designed  with  specified  risk  combinations,  the  effect  of 

chan9“  V"1  x 0,1  tUkM- 

The  second  {tart  (Section  3 ) of  this  study  involves  an  enpirical 
quantification  of  the  uncertainty  in  parameter  estimates,  and  a numerical 
determination  of  the  effect  of  such  uncertainty  on  designed  plans  and 
risks.  Since  this  requires  actual  data,  we  pursue  this  study  based  on 
9cme  of  the  data  reported  in  Schafer  (1970)  and  further  analyzed  in 
Goel  and  Joglekar  (1976  c) . 

Throughout  this  study  we  have  assumed  that  the  failure  distribution 
is  exponential  and  the  prior  for  tine  parameter  6 is  inverted  gamma. 

Only  changes  in  the  parameters  of  the  prior  distribution  have  been 
studied  in  this  report.  It  should  also  be  pointed  out  that  an  analytical 
approach  to  sensitivity  study  did  not  sen  possible  in  this  case  and  we 
had  to  resort  to  a numerical,  investigation. 
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2.  tittfo-T  OF  SPECIFIED  CHANGES  IN  PRIOR  PARAMETERS 


The  purpose  of  this  section  is  to  examine  the  effects  of  the  \m- 
certainties  associated  with  the  prior  men  and  parameter  X on  the 
designed  plan  (T,r*)  and  on  the  selected  risks.  Throughout  this  study, 
far  illustration  purposes  the  design  valuss  of  specified  MTBF  and 
miniitun  acceptable  MTBF  are  taken  to  be  100  and  50  respectively. 

The  effect  of  changes  in  T,r*,  upand  X on  various  risks  is  considered  in 
Section  2.1.  The  following  sensitivity  analyses  are  performed  in  Section  2.2. 
(1)  For  fixed  risks,  the  effect  of  small  changes  in  U^and  * on  the  de- 
signed values  of  T and  r*. 

(ii)  For  plans  designed  with  specified  risk  examination,  the  effect  of 
changes  in  U^and  X on  risks. 


2.1 . Effect  of  Changes  in  T,  r*,  Upand  X on  Risks 

The  classical  risks  are  not  influenced  lay  changes  in  g(0).  The  effect 
of  changes  in  T , r*  on  a and  B is  shown  in  Fig.  1.  The  effects  of  changes 
in  upand  X on  a and  B are  shown  in  Figures  2 (a)  and  2 (b)  while  those  of 
T and  r*  and  a and  B are  shown  in  Figures  3 (a)  and  3 (b).  Similar  results 
fo*  * and  8*  are  given  in  Figures  4 and  5. 

From  Fig.  1,  ms  see  that  when  r*  is  kept  constant,  an  increase  in  T fran 
100  to  200  hours  leads  to  a decrease  in  6 and  an  increase  in  a.  This  is 
obvious  since  larger  T far  constant  r*  inplies  larger  probability  of  rejection 
arr.  smaller  probability  of  acceptance  for  all  0.  Similarly,  far  constant  T, 
ar  increase  in  r*  from  2 to  4 implies  smaller  a and  larger  6. 


tT,r*) 

(100,4) 


Effect  of  Changes  in  fT,r»)  6n  Classical  o.C 


1.1 


100  150  200 

Specified  MTBF,  eQ 

Pig.  3(a)  Effect  of  Changes  in  T and  r*  on  a 


<T,r*) 

(100.4 


Fig.  3(b)  Effect  of  Changes  in  T and  r*  on  ? 


4(a)  Effect  of  Changes  in  u and  A on  a* 


(75,4) 


Pig.  4(b)  Effect  of  Changes  in  p 


(C  * XX  'Si 


Fig.  5(a)  Effect  of  Changes  in  T and  r*  on 


<T,r*) 


Minirnun  Acceptable  MTBF,  0 
Pig.  5(b)  Effect  of  Changes  in 


2.2  Sensitivity  of  the  Designed  Plana  and  Conputed  Risks  to  Changes  in  Prior 

Parameters. 

We  now  numerically  evaluate  the  effect  of  changes  in  u and  A on  the 

b* 

designed  T and  r*,  and  the  risks. 

2.2.1.  The  Design  Regions 

The  design  regions  being  considered  for  this  ftudy  are  as  follows.  The 
primary  design  region  considered  in  this  study  has  four  base  points  consisting 
of  all  combinations  of  Up=  75  , 125  and  A * 3,5.  Around  each  base  point, 
a factorial  design  is  constructed  with  a 5%  change  in  Upani  1 unit  change 
in  A on  either  side  of  the  base  point.  This  implies  that  in  the  subeequent 
analysis  the  effect  of  a 10%  change  in  w^and  2 unit  change  in  A will  be 
evaluated.  Additional  factorial  regions  considered  are:  (Dy^  - 120, 

140;  A « 2.5,  3.5.  (2)  yp  - 120,  135;  A » 2.5,  3.5.  (3)  « 200, 

300;  A - 2.5,  3.5  and  (4)  yp  - 187.5,  262.5;  A - 2.5,  3.5. 

It  should  be  noted  that  these  regions  have  been  rhnoni  for 
illustration  pirpoees.  Similar  analyses  could  be  dons  for  other  regions 
of  interest. 
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The  choice  of  the  changes  in  Up  and  A is  governed  by  the  uncertainties 

associated  with  the  prior  parameters.  Figures  6(a)  to  6(h)  show  the  plots 

of  the  inverted  gamma  prior  density  for  each  of  the  five  points  of  the 

eight  factorial  regions  in  this  study.  The  plots  indicate  that  the  selection 

of  changes  in  Mp  and  A is  reasonable  {except  for  designs  7 and  8 in 

Figures  6(g)  and  6(h)).  In  other  words,  if  the  prior  is  based  vpan  a 

reasonable  amount  of  data,  the  changes  in  yp  and  A considered  in  this  study 

would  be  of  the  order  of  2g*  and  2o*  respectively.  It  is  suggested  that 

* 

this  be  further  verified  by  conducting  an  analysis  of  the  real  data. 

Table  1 gives  the  values  of  the  probabilities  P(0  £ 0^,  P(e1  < 6 < 0Q) 
and  P(0  > 0Q)  for  the  forty  prior  distributions  considered  above.  Fran 
this  table  we  note  that,  in  general,  the  probabilities  in  the  tail  regions 
0 * e0  and  e - change  considerably  with  changes  in  Up  and  X.  This  in- 
dicates that  the  designed  plans  , especially  those  that  use  these  areas 
as  risk  criteria,  may  be  sensitivie  to  changes  in  the  prior  parameters. 


2.2.2 


?s  in  Designed  I Ians  For  FixSd  Risks 


For  purposes  of  this  analysis,  plans  were  designed  for  specified  risks 
and  parameters  up  and  A.  The  results  of  this  analysis  are  shown  in  Tables 
2 through  7.  In  these  tables,  the  center  point  shows  the  designed  plans 
(T,r*)for  the  specified  risk  values  and  the  corresponding  prior  parameters. 
The  values  at  the  comers  of  the  various  squares  are  the  designed  plans  if 


6(d)  inverted  Ganna  Density 
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Table  1 


Effect  of  Changes 

in  u and  x 
P 

on  Prior  Probabilities 

(6^100,6^50) 

UP 

X 

P^ep 

P(e1<e<e0) 

P(6>fc0J 

?1.25 

2 

.582 

.247 

.171 

71.25 

4 

.381 

.433 

.186 

78.75 

2 

.532 

.269 

.199 

78.75 

4 

.306 

.460 

.234 

75. 

3 

.422 

.371 

.207 

118.75 

2 

.313 

.336 

.551 

118.75 

4 

.076 

.415 

.508 

131.25 

2 

.262 

.351 

.839 

131.25 

4 

.047 

.367 

.586 

125. 

3 

.125 

.392 

.483 

71.25 

4 

.381 

.433 

.186 

71.25 

6 

.286 

.543 

.172 

78.75 

4 

.306 

.460 

.234 

78.75 

6 

.204 

.565 

.231 

75. 

5 

.285 

.509 

.206 

Table  1 (Continued) 


yp 

\ 

Pte<e1) 

P(e1<e<e0) 

P(e>eQ) 

120. 

2.5 

.206 

.380 

.414 

120. 

3.5 

.101 

.609 

.489 

140. 

2.5 

.136 

.361 

.503 

140. 

3.5 

.052 

.347 

.601 

150. 

3. 

.109 

.382 

.509 

120. 

2.5 

.206 

.380 

.414 

120. 

3.5 

.101 

.409 

.489 

135. 

2.5 

.151 

.367 

.482 

155. 

3.5 

.062 

.364 

.574 

127. 

3. 

.117 

.387 

.496 

200. 

2.5 

.035 

.247 

.717 

200. 

3.5 

.005 

.160 

.832 

300. 

2.5 

.092 

.905 

500. 

3.5 

.029 

.971 

250. 

3. 

.003 

.105 

.892 

187.5 

2.S 

.047 

.273 

.680 

»Ti 

1 

O0 
«— 4 

3.5 

.009 

.193 

.798 

202.5 

2.5 

.007 

.138 

.855 

2o2 . 5 

3.5 

.000 

.057 

.943 

225. 

3. 

.007 

.147 

.846 

25 
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Effect  of  Changes  in  w and  X on  T and  r*:  Criterion  (a,B*) 


Effect  of  Changes  in  u and  X on  T and  r*:  Criterion  (a,7f) 


Prior  Mean 


Effect  of  Changes  in  u and  X on  T and  r*:  Criterion  (a*, 8*) 
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Table  6 

Effect  of  Changes  in  u and  X on  T and  r*:  Criterion  (P(A),6*) 


(604,6) 


(387,3) 


m 

& 


P(A)-.8 

B-.05 

• 

(468,4) 


(986,12)  (456,4) 

200  250  300 


2.5 


changed  to 


the  risk  value*  were  kept  the  earn  while  the  prior  parameters  ware 
the  values  corresponding  to  the  oomer  being  considered.  Scree  of  the 
observations  from  these  Tables  are  sxxmarized  below.  It  should  be  noted  that 
the  inferences  drawn  from  these  tables  are  dependent  on  the  design  region 
being  investigated.  Similar  analyses  should  be  performed  for  regions  that 
may  be  of  interest  for  a different  problem. 

Fixed  a and  B*  (Table  2) 

Factorial  I 

As  }ip  increases,  T decreases.  The  decrease  is  larger  for  higher  X; 
r*  is  almost  uiaffectad. 

As  X increases,  T decreases.  Ihe  decrease  is  larger  for  higher 
r*  is  almost  unaffected. 

As  tip  end  x increase  simultaneously,  T decreases. 

Factorials  II,  III,  IV 

As  the  prior  beoomes  favorable,  test  times  reduce  and  in  four  cases, 
via.  4ip#X)  - (118.75,4),  (131.25,4),  (118.75,6)  and  (131.25,6),  acceptance 
is  permitted  without  any  tasting  for  values  of  a and  (F  shown  in  the  Table. 
Fixed  5 and  f(Table  3) 

Within  the  region  considered,  an  increase  in  Up  iaplies  an  increase  in 
T and  r*.  An  increase  in  X also  leads  to  an  increase  in  T and  r*. 
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Fixed  g*  and  8*  (Table  4) 

Depending  upon  the  values  of  yp  and  X,  and  increase  in  Up  may  lead  to 
an  increase  or  a decrease  in  T and  r*.  Similarly,  an  increase  in  * may 
lead  to  larger  or  smaller  T and  r*. 

Fixed  5 and  0*  (Table  5) 

An  increase  in  w leads  to  a reduction  in  T and  r*  for  constant  A. 

P 

An  increase  in  A may  yield  larger  or  smaller  T and  r*  depending  upon  the 
parameter  values. 

Fixed  P(A)  and  g*  (Table  6} 

An  increase  in  yp  and/or  A leads  to  a reduction  in  T and  r*. 

Fixed  P (A)  and  6 (Table  7) 

An  increase  in  yp  and/or  A reduces  T and  r*. 

2.2.3  Changes  in  Ccrputed  Risks  For  Fixed  Plans 

The  design  region  described  in  Section  2.2.1  was  considered  for 

numerically  studying  the  effects  of  changes  in  yp  and  A on  the  risks. 

The  risks  were  evaluated  at  the  four  corner  points  of  the  design  assuring  that 

the  plan  designed  for  the  center  point  is  also  enployed  at  the  corresponding 

comer  points.  The  results  are  shown  in  Tables  8 to  13  and  the  effects  of 

u and  A on  the  risks  can  be  observed  from  rinses  tables. 

P 


Effect  of  Changes  in  u and  > on  Risks:  Criterion  (a#B+) 


Prior  Mean 


Prior  Mean 


Effect  of  Changes  in  u and  X on  Risks:  Criterion  (a*«6*) 


Prior  Mean 


Effect  of  Changes  in  u and  X on  Risks:  Criterion  (a,6#) 


, Effect  of  Changes  in  y and  X on  Risks:  Criterion  P(R),6* 


Prior  Mean 


2.2.4  Sane  General  Cormanta 


In  the  following  paragraphs,  we  provide  sane  brief  ocranents  based  on 
the  analyses  in  Sections  2.2.2  and  2.2.3.  As  was  pointed  out  earlier, 
such  inferences  are  true  for  the  desi^i  region  considered  and  should 
not  be  construed  to  be  true  in  general. 

(a)  All  the  plans  considered  are  quite  sensitive  to  changes  in  the  prior 
panmeter  values,  both  in  terns  of  T and  r*  and  in  terms  of  the  risks. 
The  CP  (A) , 8*)  criterion  shows  the  largest  sensitivity.  For  example, 
in  Table  6 we  see  that  foryp  ■ 120  and  A - 2.5  we  have  T = 1134  and 
r*  ■ 23.  For  Pp  ■ 135  and  A ■ 3.5,  the  design  values  drop  to  T = 16 
and  r*  « 0. 

(b)  While  evaluating  the  effects  of  changes  in  the  prior,  the  changes  in 
T and  r*  must  be  jointly  considered  since  a part  of  the  change  in  T 

is  due  to  the  change  in  r*.  For  exanple,  for  the  case  just  considered, 
if  r*  is  changed  frcra  0 to  23,  T will  be  considerably  larger  than  16 
and  will  approach  1134. 

(c)  For  a proper  perspective  on  the  sensitivity  analysis,  we  must  look  at 
the  changes  in  T,  r*  and  risks  jointly.  Using  the  (P  (A) , p*)  criterion, 
let  us  consider  an  example  to  explain  this  point.  Let  us  assume  that 
the  true  prior  has  Mp  ■ 120  and, X ® 2.5  and  we  wrongly  assume  it  to  have 
^p  * 127.5  and  A-  3.0.  For  specified  P(A)  « 0.8,  8*  « 0.045  and 

9q  * 100,  from  Table  6 we  have: 
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Assumed  Prior:  Designed  Plan  (238,  4) 
True  Priori  Designed  Plan  (1134,  23) 


A. 


Thus,  the  error  in  the  prior  has  led  to  reduced  T and  r*.  Once  the 
plan  (238,  4)  is  implemented,  the  true  risks  will  correspond  to  the  true 
prior  and  will  not  equal  0.8  and  0.045. 

From  Table  12,  we  have  the  corresponding  risks  as: 

Assumed  Prior:  P(A)  = 0.8,  8*  * 0.045 
True  Risks:  P(A)  * 0.7303,  8*  = 0.0752 

This  example  shows  that  the  error  an  the  prior  has  led  to  a reduced  T and  r* 
and  higher  risks.  Similarly  it  may  happen  that  T and  r*  get  larger  and  the 
true  risks  get  smaller. 

It, therefore, follows  that  if  we  consider  the  cost  of  testing  and  the 
oost  of  making  wrong  decisions,  an  error  in  the  prior  may  lead  to  increased 
testing  cost  and  reduced  oost  due  to  risks  or  vice  versa.  This  balancing 
of  costs  indicates  that  the  oost  may  be  less  sensitive  to  changes  in  the 
prior  and  points  out  the  need  to  consider  the  oost  structure  in  the  design 
of  plans. 

2.3  Teat  Times  and  Test  Tine  Relationships 

The  effect  of  varying  ^ and  varying  X on  T is  shoU^an  Figures  7 and  8 
respectively  for  several  different  criteria.  The  slopes  of  these  curves 
are  a measure  of  the  sensitivity  of  T to  changes  in  ^ and  X.  A 
comparison  of  test  times  based  an  same  nunerlcal  values  of  risks  is  somewhat 
meaningless  since  the  design  criteria  are  different.  The  plots  in  Figures 
7 and  8 are  provided  only  for  the  sake  of  acnpletenass. 
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NOTE: 


O 2.0 

Prior  Mean  vp  X eQ 

Figure  7:  Effect  of  ^ on  T»  (K  ■ 2,X  ■ 4) 

Discontinuities  exist  on  all  curves;  curves  are  drawn  only  for 
information.  A9 
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3.  QUANTIFICATION  CF  UNCERTAINTY  IN  PRIOR 
PARAMETER  ESTIMATES 

In  the  previous  Section  we  studied  the  sensitivity  of  designed  plans 
and  ocnputed  risks  to  specified  changes  in  the  prior  distribution  in- 
dependent of  the  availability  of  data.  Near  we  consider  the  situation 
whan  the  prior  parameters  are  obtained  from  failure  data. 

Since  the  parameters  are  estimated  fron  failure  data,  the  uncertainty 
in  estimation  can  be  quantified  by  the  varianoe-oovariance  matrix  of 
X,  Y as  described  in  Goel  and  Joglekar  (1976  c) . If  the  correlation 
coefficient  p<J^a  between  X and  9/  is  zero,  the  sensitivity  studies  can  be 
carried  out  by  trusting  3 and  y to  be  independent.  However,  if 

9\  * + 0,  the  joint  aonfidanoa  regions  for  X and  y are  needed.  Since  these 
A»Y 

regions  for  —11  amples  cannot  be  obtained  mathematically  for  the  present 
study,  wa  resort  to  altamativa  approaches. 

The  first  approach  is  to  obtain  such  regions  by  simulation  and  to  use 
the  points  along  the  100  (l-**)  I confidence  contour  to  provide  a measure 
of  uncertainty  in  the  parameter  estimation. 

The  seoond  approach  is  to  obtain  the  likelihood  plots  of  X and  y based 
on  observed  data  and  use  them  to  obtain  a measure  of  uncertainty.  The 


third  approach  is  too  study  the  asynpbotic  bivariate  normal  distribution 

-A  A 

of  X,  and  y. 


~ '.'it  -l  . 
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In  the  following  subjections  we  use  these  three  approaches  for 
sensitivity  studies.  TWO  sensitivity  analyses  are  conducted.  First  the 
effect  of  the  uncertainty  in  the  prior  parameter  estimates  on  the  normalized 
test  time  T*  («  T/8Q)  and  aooaptanoe  number  r*  for  a single  sanple  truncated 
plan  for  a system  Is  determined.  Four  risk  criteria,  namely  (a ,8), 

(a, 8*),  (a*,0*)  and  (P(R),$*)  (see  Goel  and  Joglekar  1976(b) far  details) 
are  considered.  In  each  case  designed  (T*,r*)  valves  are  obtained  for 
specified  risk  values  and  estimated  (\,y) . Then,  keeping  the  risks  constant, 
designed  values  are  obtained  for  same  selected  values  of  \ and  y using  one 
of  the  above  mentioned  three  approaches.  The  variation  in  the  designed 
(T*,r*)  values  provides  a measure  of  sensitivity. 

A second  sensitivity  analysis  is  conduct*!  to  determine  the  effect  of 

A A 

uncertainty  associated  with  1 and  y on  the  risks.  The  four  risk  criteria 

given  above  are  considered.  In  each  case,  designed  T*  and  r*  values 
\ A 

corresponding  to  *and  Y are  used  to  ocnpute  appropriate  risks  at  selected 
(*  ,Y ) points.  The  variation  In  risks  constitutes  a measure  of  sensitivity. 

The  reasons  behind  the  observed  sensitivity  and  the  choice  of  criteria 
an  the  basis  of  sensitivity  are  discussed  next.  An  approach  to  determine  the 
amount  of  data  necessary  to  obtain  a prespecified  sensitivity  is  outlined. 


Also,  a oonparison  of  the  results  from  the  simulation,  likelihood  and  bi- 
variate normal  approaches  provides  practical  guidelines  for  the  quantification 
of  uncertainty. 


f 


3.1  Variance  Covariance  Matrix  for  X ,y 

Parameter  estimation  has  been  considered  in  detail  in  Goel  and  Joglekar 
(1976  c) . It  has  been  shown  that  if 

f (t| © ) - i e-t^  ; ti  0,  6 >0  (3) 

and 


0* 

A 


g(0)  «-* — 0 U X)  en/B  ; 0,y,A  > 0 (4) 

r cx> 

then  the  marginal  density  of  the  mxtber  of  failures  X in  a fixed  time  T 
is  negative  binomial  i.e. 

*<x>  • ITTPlit  x - O.I.2.3.,.  (5) 

The  convene  is  also  true.  If  the  observed  frequency  nx  corresponding 

to  x - 0 is  not  reaordsd  ben  X has  a truncated  negative  binomial  dis- 
tribution given  by 


f (x)  - -JE+gil* . (^x  {£-))  x - 1,2..  16) 

xitt-DUl-^*) 

It  should  be  noted  that  the  prior  parameters  X and  y also  form  the 
of  the  distribution  of  the  aba  enable  random  variable  X.  The 
parameters  can  be  estimated  fay  fitting  a negative  binomial  or  a truncated 
negative  binomial  distribution  to  the  observed  values  of  X dtpending  on 
whether  observations  corresponding  to  x ■ 0 are  recorded  or  not.  In  our 
case  observations  for  x - 0 are  missing  and  hence  the  maximun  likelihood 


*nL  ■ i n^UntA+x-l) ! -In  (x) ! -fcn(X-l) ! + in 


Estimates  x and  y can  be  obtained  by  an  iterative  solution  of  inl/BX  - 0 
and  InL/frf  ■ 0.  The  variance  oovarianoe  matrix  for  X and  y can  be 
is— ri rally  obtained  by  taking  the  seoond  partial  derivatives  of  inL  and 
is  given  by. 


covCY,*)  var(X) 
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tasting  and  is  the  observed  frequency.  Since  observations  corresponding 
to  x • 0 are  missing,  a truncated  nsgative  binomial  is  fitted  and  the 
maxlntm  likelihood  estimates,  their  variances  and  the  correlation  coefficient 
are: 


System 

A 

X 

A 

y 

A 

A 

0* 

Y 

Ky 

Tranmaitter 

0.2054 

680.0 

0.28 

318.0 

0.86 

Search  Indicator 

0.4407 

1152.0 

0.40 

542.0 

0.90 

Search  IWP8 

0.8646 

714.0 

0.31 

2U.0 

0.89 

Since  is  approximately  0.9  for  all  three  aaaes,  the  individual  variance 
eetimetee  do  not  provide  sufficient  quantification  of  uncertainty  and  joint 
confidence  regions  for  A and  y bus*  be  obtained.  As  mantionsd  earlier, 
explicit  expression  for  the  joint  density  does  not  seem  feasible  since  the 
parameter  eetimetee  themselves  are  not  explicitly  known. 


3.2  Simulation  Study 

A study  was  conducted  to  obtain  the  joint  and  the  marginal 

densities  of  X and  y.  The  study  can  best  be  explained  by  an  example. 
Consider  the  failure  data  in  Table  1 for  the  search  indicator.  Based  upon 
55  failures  the  estimated  parameters  are  *X  « 0.4407  and  y • 1152.0. 

The  steps  in  the  simulation  study  are  as  follows: 

1.  Take  X ■ 0.4407  and  y - 1152.0  and  generate  55  obeervations  from  a 
truncated  negative  binomial  distribution  with  parameters  \,y. 

2.  Estimate  the  parameters  X and  y for  the  generated  failure  data  using 

likelihood  estimation  procedure. 

A /\ 

3.  Repeat  the  above  steps  N (say  200)  times.  A plot  of  X vs  y gives  the 

A A 

joint  density.  The  histograms  for  X and  y constitute  the 
simulated  marginal  densities. 

i 

Using  the  above  procedure,  the  following  four  simulationsswere 
oorauctea* 


System 

No.  of 
Observations 
for  each  simulation 

No.  of 
Simulations 

Transmitter 

0.2054 

680.0 

50 

200 

Search  Indicator 

0.4407 

1152.0 

55 

200 

Search  WPS 

0.8646 

714.0 

74 

200 

Search  Indicator 

0.4407 

1152.0 

55 

1000 

52 
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Based  on  the  above  simulations,  four  sets  of  X and  y values  are 
obtained,  the  first  three  consisting  of  200  pairs  of  values  and  the  fourth 
one  with  1000  pairs  of  X,y.  The  first  three  sets  of  X,y  values  are  given 
in  Appendix  A,  Tables  A.l  to  A.  3,  The  simulated  joint  densities  and 
the  simulated  marginal  densities  are  obtained  from  these  sets  of  (X,y) 
values.  The  joint  densities  for  cases  (a)  through  (c)  are  given  in 
Figures  9 to  11  respectively,  and  the  marginal  densities  for  X and  'y  for 
the  four  cases  are  plotted  in  Figures  12  to  19.  The  cmpited  values  of 
sanple  mean,  sanple  variance,  skewness  0^  and  kurbosis  are  *i<r>  given 
in  Figures  12  to  19. 

/>  * 

The  joint  density  of  X and  y has  the  interpretation  that  the  95% 
contour  corresponds  to  the  95%  confidence  region  for  X and  y.  Points 
along  the  95%  confidence  contour  can  be  judiciously  selected  to  quantify 
the  uncertainty  associated  with  the  prior  panmeters  for  the  purpose  of 
sensitivity  analysis. 

The  following  observations  may  be  made  f ran  the  above  simulation  results, 
(i)  Sens  of  the  estimated  parameter  values  are  negative  and  are  in- 
atfcviasible.  These  negative  values  may  be  caused  because  of  the 
nature  of  the  simulated'  data  in  these  few  cases.  This  is  not 
•uprising  when  an  iterative  procedure  is  used  to  get  the  maxinun 
likelihood  estimates,  especially  when  the  value  of  X is  so  close  to 
*«eo.  <*«  ocrrncnly  used  approach  is  to  take  the  negative  value  to  be 
•lightly  greater  than  aero.  Another  possibility  is  to  use  Brass's 
modified  manant  estimates.  These  approximations  are 
only  because  it  is  known  that  the  data  is  fran  a TKBD. 
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CliM  Interval 

rifMra  14.  Marginal  Omalty  of  i(».  0.4407,  t-  1157.0) 


Class  Interval 

Marginal  Density  of  y(Y-0.M07,  y-1152.0) 
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Pt(ur«  16.  Marginal  Ocnalty  of  i(i>0.H626,  v-714.0) 


Class  Interval 

Marginal  Density  of  y(X-0.«626.  v-714.0) 
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(11)  The  aanputed  averages  X and  y are  aonslstantly  larger  than  X and  y 
indicating  the  snail  sanple  bias  of  naxiitun  likelihood  estimates. 

The  bias  remains  the  same  for  200  simulations  and  1000  simulations 
and  depends  upon  the  nimfcer  of  observations  in  each  simulation 
and  not  an  the  mmfeer  of  simulations. 

(iii)  The  small  sanple  average  standard  deviations  $£  and  are  larger 
than  the  standard  deviations  obtained  from  the  information  matrix. 

(iv)  The  marginal  distributions  of  X and  y are  skewed  to  the  ri$it  as 
also  indicated  by  the  positive  values  of  skewness. 

(v)  The  marginal  distributions  are  generally  leptbkurtic.  200 
simulations  are  insufficient  to  estimate  as  evidenced  by  the  fact 
that  (%  * 13.6  for  200  simulations  and  a^(X)  « 3.37  for  1000 
sinulations,  for  cases  (b)  and  (d).  Figures  14  and  18  respectively. 
This  is  expected  since  estimation  of  higher  moments  will  require 
larger  simulation. 

(vi)  200  sinulations  appear  to  be  sufficient  to  estimate  the  joint  and 
the  marginal  densities  since  these  densities  remain  essentially  the 
seme  for  1000  simulations. 

Results  of  sensitivity  analyses  are  discussed  later  in  Section  3.5. 
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Figure  22.  Contours  of  Normalized  Likelihood  (Search  MVPS  Data,  \ = 0.8646,  y = 714) 


I 


The  likelihood  contours  do  not  have  the  confidence  region  interpretation. 
However,  an  empirical  comparison  of  the  likelihood  contours  and  the 

empirical  joint  densities  shows  that  the  5%  likelihood  contour  reasonably 

% 

approximates  the  95%  confidence  region  and  may  be  taken  to  be  an  easy  to 
use  approximate  quantification  of  uncertainty. 

Sensitivity  analyses  based  on  these  contours  are  discussed  later  in 
Section  3.5. 


3.4  Asyaptotic  Bivariate  Normal  Contours 

Maximum  likelihood  estimates  ace  asymptotically  normally  distributed. 

A A 

Therefore,  a and  Y have  an  asynptotic  joint  bivariate  normal  distribution 
given  by 


”\)2 


x,r  la -r 

-2Pt  „£*>  (F»  + <jp>2» 
X'y 


*,y 


) * 


(12) 


x y y 

Centaurs  containing  100  (1-a)  % of  the  distribution  are  given  by 

c|r>2  • 2 ox,;  &r>  & * &2 ' ‘2(1'°;  ;>  l0‘(1-a)  (13) 

X X y y ,T 


Figures  23,24  and  25  show  plots  of  such  contours  f or  a * 0.05,  010  and 
0.1587  for  the  three  cases  being  considered.  The  plots  are  drawn  by 
taking  0^,0*  and  to  be  numerically  equal  to  a^,o^  and  a 
respectively. 

These  contours  can  be  interpreted  as  asymptotic  confidence  regions 
for  X and  y.  A oonpariaon  with  the  onall  sanple  results  by 

siaulation  indicates  that  these  contours  may  be  considered  to  form  a 
quantification  of  uncertainty. 
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Contours  of  Asymptotic  Normal  Distribution  of  ( A , Y) 
(Search  Indicator  Data) 


The  three  approaches  discussed  above  were  intended  to  obtain  same 
quantification  of  the  prior  parameters  X and  y.  Using  the  results  from  these 
approaches  of  analysis,  the  following  two  types  of  sensitivity  were 
conducted. 

1.  For  each  of  the  three  examples  under  consideration,  X and  y values 
along  the  95%  confidence  contour  were  judiciously  selected.  For  the 
four  risk  criteria,  namely,  (a*, 6*),  (a,B),  (a, 6*)  and  (P (R) , B*) , 
plans  were  designed  to  obtain  T*  and  r*  at  the  selected  parameter  values 
for  two  discrimination  ratios,  K = 2 and  K * 1.5.  The  results  are 
given  in  Tables  15  to  20. 

A ^ 

2.  For  each  of  the  four  risk  criteria,  plans  were  designed  at  X and  Y. 

Using  the  designed  T*  and  r*  values,  the  risks  were  computed  at  each 
of  the  X,y  values  selected  above.  The  "results  from  this  analysis  are 
tabulated  m Tables  21  to  26. 

The  first  sensitivity  analysis  gives  the  effect  of  the  uncertainty 
associated  with  the  parameters  on  the  designed  plans  and  the  second  sensitivity 
analysis  gives  the  effect  on  risks. 

The  following  observations  are  made  from  the  results  of  the  above 
sensitivity  analyses  as  given  in  Tables  15  to  26. 

(l)  In  general , tne  designed  plans  and  risks  are  significantly  sensitive 
to  changes  an  the  prior  parameters  for  the  data  under  consideration. 

$ Changes  in  T*  by  a factor  of  4,  changes  in  r*  by  a factor  of  2 
and  changes  in  risks  by  a factor  of  3 are  faitly  common. 
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TABLE  A.1 


SIMULATED  VALUES  OF  1 AND  Y (*-0.2054.  Y-680) 


GAMfA 

LAMBDA 

GAMMA 

LAMBDA 

1626.2998 

1.2906 

470.2000 

-0.0349 

1163.7000 

0.6500 

454.8999 

0.1662 

1188.8999 

0.5156 

800.3999 

0.2258 

338.8999 

-0.0752 

887.0000 

0.4657 

1053.3999 

0.4376 

621. 7998 

0.1914 

644.5999 

0.1307 

1199.3999 

0.5029 

488.2000 

-0.0717 

671.3999 

0.1741 

1404.5999 

0.5922 

261.3999 

-0.2296 

819.0000 

0.2997 

643.7998 

0.1623 

885.2000 

0.4050 

507. 7998 

0.0958 

1298.2998 

0.6401 

587.7000 

0.1637 

412.7998 

-0.1883 

594.0000 

0.1707 

978.3999 

0.7273 

1430.0000 

0.5358 

1383.0000 

0.6035 

386.2998 

-0.0608 

730.5000 

0.4183 

1308.7998 

0.9830 

906.8999 

0.1529 

1094.8999 

0.2526 

691.5000 

0.2390 

394.2000 

-0.1483 

555.0000 

0.2494 

817.2000 

0.2079 

1254.5999 

0.6296 

1160.3999 

0.3279 

785.2998 

0.1398 

810.7000 

0.1061 

723.2000 

0.2949 

573.5000 

0.0935 

748.2000 

0.2448 

2302.5000 

1.5078 

604.2000 

0.1214 

374.7998 

-0.1243 

655.5000 

0.2488 

1044.8999 

0.9630 

793.0999 

0.2512 

1107.5999 

0.5892 

769.3999 

-0.0403 

905.3999 

0.0970 

734.0000 

0.1953 

566.2998 

0.2369 

383.0000 

-0.1965 

1267.0000 

0.3970 

713.7998 

0.0072 

451.5000 

0.0656 

1849.5000 

1.5158 

687.2000 

0.2565 

916.5000 

0.6314 

522.8999 

0.0355 

809.5999 

0.3364 

853.7000 

0.4314 

660.3999 

0.0665 

1798.3999 

0.6599 

763.0000 

0.3659 

589.7998 

0.1855 

2387.8999 

1.6382 

499.5999 

0.1583 

532.7000 

-0.0325 

2011.3999 

0.8416 

762.0000 

0.4822 

326.5999 

0.0645 

798.7998 

0.3744 

719.0999 

0.4675 

936.0000 

0.5954 

259.0999 

-0.3032 

5323.8984 

2.6631 

447.7000 

0.0369 

754.8999 

0.1902 

1302.2998 

0.7074 

1117.0999 

0.5014 

705.0999 

0.3410 

684.0999 

0.6191 

754.2998 

-0.0156 

719.0000 

0.1013 

819.2998 

0.3208 

545.0999 

0.1531 

700.7998 

0.2490 

415.5000 

0.0281 

1563.0000 

1.0285 

480.8999 

-0.1029 

487.2000 

-0.0259 
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TABLE  A.l(Contlnu«d) 


GAMtA 

LAMBDA 

GAMMA 

LAMBDA 

1556.7000 

0.8256 

692.0000 

0.0150 

1962.2998 

1.0254 

529.7998 

0.1260 

809.5000 

0.2210 

592.7998 

0.0888 

640.0000 

0.2772 

646.8999 

0.2082 

611.0999 

0.0703 

1071.5000 

1.2839 

1084.5000 

0.9083 

821.8999 

0.2456 

1111.2000 

0.4549 

1266.7998 

0.5211 

992.2000 

0.3611 

567.0000 

0.0370 

1031.5000 

0.4930 

1153.7998 

0.4305 

510.7000 

0.0270 

694.5999 

-0.0950 

3705.2000 

1.9025 

395.0999 

0.1135 

474.7000 

0.0098 

962.5999 

0.2435 

1242.5000 

0.7404 

860.0999 

0.4695 

238.6000 

-0.2872 

816.2000 

0.0748 

971. 7000 

0.3659 

685.0999 

0.3027 

970.7000 

0.2967 

569.0999 

-0.0184 

1855.0000 

0.9946 

1208.3999 

0.4195 

331.8999 

-0.2345 

423.7998 

-0.0040 

827.0999 

0.5062 

360.2000 

-0.0028 

697.0000 

0.1770 

557.2998 

0.0563 

765.0999 

0.2246 

514.5000 

-0.0260 

385.0999 

-0.1986 

607.2998 

0.2537 

1305.5000 

0.2152 

753.5999 

0.1575 

743.7000 

0.3842 

494.3999 

0.0407 

1376.5999 

0.4443 

615.0999 

-0.0560 

582.5999 

0.4101 

658.2000 

0.3234 

464.0000 

0.0483 

768.0999 

0.1770 

2172.5000 

1.0529 

475.7000 

0.0915 

920.5999 

0.3134 

516.5000 

0.1115 

472.7998 

0.0673 

802.7000 

0.3845 

555.0999 

0.1737 

765.8999 

0.1296 

346.3999 

-0.2799 

1188.2998 

0.4222 

848.7000 

0.2901 

433.0000 

0.2708 

534.7000 

-0.0929 

452.0000 

-0.0588 

436.0999 

0.0148 

1201.7998 

0.6266 

1007.8999 

0.5750 

2455.5999 

1.5037 

302.7000 

-0.0306 

1207.2000 

0.6057 

608.5999 

0.2744 

717.2000 

0.2996 

843.5000 

0.4082 

401.7998 

-0.0194 

595.7998 

0.2308 

838.7998 

0.1926 

943.7000 

0.3042 

810.5000 

0.2416 
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TABLE  A.l(Contlnu«d) 


GAWA  LAMBDA 

948.2998  0.0787 

808.8999  0.4974 

887.8999  0.3338 

2238.0999  1.7486 

608.7998  0.2328 

729.7998  0.4858 

1030.2998  0.5279 

534.0999  0.2690 

553.3999  0.1112 

527.8999  0.1504 

961.7000  0.3197 

771.5999  0.1280 

915.2998  0.1516 

439.7998  -0.1480 

505.0999  0.0980 

950.0000  0.5027 

1414.5000  0.3388 

530.8999  0.1046 

415.7000  0.0260 

989.5999  0.3435 

438.3999  0.1724 

1298.0999  0.5174 

327.0999  -0.1732 

882.8999  0.635/ 


TABLE  A. 2 


SIMULATED 

VALUES  OF 

A AND  Y (A-0.4407.  v-1152) 

GAMMA 

LAMBDA 

GAMMA 

LAMBDA 

1702.2000 

0.6841 

737.0000 

0.1345 

467.7000 

0.0101 

1567.2000 

0.3940 

780.2998 

0.4946 

1081.5999 

0.3940 

2371.3999 

1.1299 

1379.0000 

0.6238 

1710.2998 

1.0923 

639.7998 

0.2651 

1303.7998 

0.8202 

789.2998 

0.4063 

900.5000 

0.2266 

886.2000 

0.4401 

561.5999 

0.0343 

1517.8999 

0.8203 

435.3999 

-0.2675 

1287.3999 

0.4815 

431.3999 

-0.1738 

650.7000 

0.1680 

1016.0000 

0.3795 

1457.2998 

0.3480 

895.2998 

0.0839 

1937.7000 

1.0283 

1004.0000 

0.2978 

1426.0999 

0.3494 

1599.3999 

0.5365 

1539.0999 

0.4427 

1434.7000 

0.5740 

3451.8999 

2.0869 

918.0999 

0.3855 

1398.0999 

0.6663 

1166.5999 

0.4846 

828.0000 

0.2882 

820.2000 

0.2692 

1375.5999 

0.8634 

1423.0000 

0.4880 

1449.7000 

0.8176 

3209.8999 

1.6547 

845.5999 

0.2560 

1114.8999 

0.5752 

409.7998 

-0.1160 

1358.0999 

0.7036 

1012.0000 

0.0728 

1256.0000 

0.5999 

647.5000 

-0.0392 

1129.7998 

0.5222 

873.5000 

0.1372 

2508.0000 

1.2068 

948.5999 

0.4588 

828.8999 

0.3419 

1051.7998 

0.3737 

1580.0999 

0.8331 

1127.5999 

0.4132 

1136.5000 

0.4125 

3291.2998 

2.0623 

1140.5999 

0.2942 

1310.2000 

0.6711 

836.7000 

0.1015 

874.2998 

0.4280 

1004.7998 

0.3558 

5752.0977 

3.2755 

1917.0000 

0.4030 

1151.2998 

0.3762 

1093.2998 

0.2351 

1648.0999 

1.2471 

1455.0000 

0.4605 

1501.2998 

0.6825 

941.5000 

0.2675 

1001.2998 

0.1454 

2580.3999 

1.2396 

1288.2998 

0.3462 

1475.5999 

0.4539 

1431.0000 

0.6838 

2319.5000 

1.4120 

1207.5000 

0.3858 

604.0999 

-0.0280 

913.2998 

0.1754 

1096.5999 

0.5077 

1567.3999 

0.6419 

2516.0000 

1.4802 

1842.5000 

1.1328 

650.5000 

0.0346 

550.5000 

-0.0712 

1107.8999 

0.2620 

1424.8999 

0.6789 

1710.5999 

0.8435 

1064.8999 

0.1659 

1556.2000 

0.7049 

2475.5000 

1. 2642 

2113.5000 

1.0854 

1035. 7998 

0.5156 

1139.3999 

0.3831 

2233.0000 

0.8375 

4418.0000 

2.2437 

569.5999 

0.0236 
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TABLE  A.  2 (Continued) 


GAMMA 

LAMBDA 

GAMMA 

LAMBDA 

1071.5000 

0.3751 

2333.7000 

1.2734 

1586.8999 

0.5064 

1048.7000 

0.4293 

793.7000 

0.1489 

1069.7000 

0.3281 

1822.5000 

0.7558 

610.0999 

0.0235 

1854.5999 

0.9419 

1109.7000 

0.2634 

893.0000 

0.2793 

944.7998 

0.3248 

1399.8999 

0.8775 

759.8999 

-0.0230 

605.2000 

-0.0116 

1171.0999 

0.7804 

521.5000 

0.0624 

794.0000 

0.3408 

789.7000 

0.2805 

1021.2000 

0.4678 

2193.0999 

0.7700 

847.5000 

0.1100 

1214.3999 

0.4833 

2632.5000 

1.5680 

1852.2998 

0.8940 

1517.8999 

0.7728 

2335.2000 

0.9891 

977.0000 

0.4338 

475.0000 

0.0224 

1844.2000 

0.4054 

864.5000 

0.1368 

1084.3999 

0.5113 

1642.7000 

0.8328 

8633.7969 

4.9344 

1313.2998 

0.7155 

828.5000 

0.1598 

957.2998 

0.2597 

1599.7000 

0.7686 

2903.5999 

1.5218 

1318.3999 

0.5314 

1301. 2998 

0.5616 

2176.0000 

0.5675 

1035.0999 

0.4596 

1131.3999 

0.5977 

1867.0000 

0.6135 

1485.8999 

0.6317 

2777.2998 

1.2009 

2743.8999 

1.7140 

1176.7998 

0.5933 

1637.5000 

0.8287 

2419.0999 

0.8406 

1275.0000 

0.3644 

1899.8999 

0.9993 

1339.5999 

0.4568 

770.7998 

0.2453 

2236.7998 

1.5002 

911.2998 

0.3599 

5005.7969 

2.0796 

1120.5999 

0.6521 

1140.0999 

0.5463 

565.5999 

0.0567 

966.7000 

0.3438 

1953.0000 

0.9921 

999.7000 

0.5474 

1597.3999 

0.8965 

1766.5000 

1.1106 

1996.5999 

0.7723 

1847.7998 

1.3600 

1766.3999 

1.1518 

2491.5000 

1.6316 

1039.0000 

0.5322 

804.8999 

0.4500 

909.0999 

0.0794 

643.2000 

0.0280 

642.2998 

0.0902 

874.2000 

0.2370 

1401.3999 

0.5301 

1192.0000 

0.5453 

1719.8999 

1.0397 

801.0999 

0.3201 

1769.2998 

1.0396 

3615.7998 

2.6707 
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A 


GAttlA 

1105. 6999 

1685.8999 

857.2000 
939.7000 

724.2000 
883.0000 

1352.2998 
1245.3999 

536.3999 

1256.2998 

1597.7998 
841.2998 

1090.5000 

1083.5000 

2995.2998 

1633.7998 

1458.5000 

3508.2998 

1290.8999 
1283.0999 

1285.2998 
1011.5999 


TABLE  A.2(Contlnued) 
LAMBDA 


0.3723 

0.7843 

-0.1207 

0.2314 

0.0179 

0.5012 

0.4283 

0.5251 

-0.1083 

0.5000 

0.8075 

0.0163 

0.3906 

0.5674 

1.7051 

0.9467 

0.7615 

2.3385 

0.7603 

0.6805 

0.6579 

0.3402 


96 


TABLE  A,  3 


SIMULATED  VALUES  OF  X AND  y (1-0. 8626.  y-714) 


GtfMA 

LAMBDA 

gamma 

LAMBDA 

1129.7998 

1.2922 

664.8999 

0.7818 

698.5000 

0.6776 

1241.8999 

1.3697 

691. 7000 

0.6410 

709.2998 

1.0558 

794.0000 

1.0091 

801.5000 

0.6861 

1044.0999 

1.3207 

667.2998 

0.6105 

872.0000 

1.0300 

666.2998 

0.7072 

913.7998 

1.2115 

905.0999 

1.0029 

999.0999 

1.2795 

668.0000 

0.5642 

565.5000 

0.7024 

861.2998 

1.2500 

795.2998 

1.0477 

970.7000 

1.2543 

767.5000 

0.9830 

449.0000 

0.3235 

792.0000 

1.0922 

1087.0999 

1.0042 

1237.5999 

1.5421 

1543.5999 

1.8180 

818.7998 

1.0479 

912.5000 

1.2924 

325.0000 

0.1347 

723.0000 

0.9620 

885.0000 

0.8389 

806.0000 

0.8293 

1129.2998 

1.3587 

805.7998 

0.9857 

1652.2000 

1.5372 

1177.2998 

1.1297 

635.8999 

0.5439 

600.2998 

0.6387 

1048.0000 

1.3599 

725.2000 

0.7326 

954.2000 

1.0195 

2072.7000 

2.3643 

1249.2998 

1.6217 

977.0999 

1.0214 

699.2000 

0.9969 

699.5000 

0.7620 

669.7998 

0.8840 

859.0999 

1.0550 

1037.7000 

1.1893 

537.7000 

0.4938 

1007.3999 

1.3464 

910.2000 

1.4206 

1148.3999 

1.6338 

1707.8999 

1.8741 

963.7998 

1.1645 

813.2998 

0.9602 

722.2000 

1.0797 

821.7998 

1.0376 

809.0000 

0.9003 

687.0999 

0.8009 

581.3999 

0.6610 

705.8999 

0.9951 

792.5999 

1.3197 

783.5999 

1.0594 

699.7000 

0.7048 

697.7998 

0.9747 

806.3999 

1.2531 

1137.2000 

1.4718 

541.8999 

0.6337 

709.2000 

0.7618 

884.3999 

1.4058 

694.3999 

0.7746 

825.0999 

0.6945 

922.3999 

1.4652 

1507.5000 

1.6750 

945.5000 

1.0305 

1047.7000 

1.0557 

1157.7000 

1.3907 

858.7998 

0.9286 

805.2000 

0.8406 

864.5999 

1.5309 

897.8999 

1.1198 

790.2998 

0.9513 

1298.5999 

1.6937 

1050.7998 

0.8512 

973,7000 

1.2667 

772.2998 

0.9723 

963.2000 

1.0230 

1307.7000 

1.6148 

638 . 7000 

0.8663 

747.5999 

0.7406 

1046  . 5000 

1.2430 

544.5000 

0.4425 

702  . 0999 

0.8391 
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TABLE  A»3(Contlau«d) 


GgtfA 

LAMBDA 

GAMA 

LAMBDA 

671.099*9 

1.0229 

1131.7000 

1.2632 

763.2998 

0.8058 

785.7000 

0.9656 

747.5000 

0.7373 

664.0000 

0.8467 

635.7998 

0.6454 

790.0000 

0.9736 

1263.0999 

1.4755 

513:2000 

0.3631 

987.0999 

1.0228 

1068.2998 

1.2913 

754.0999 

0.9719 

729.2000 

0.6888 

1119.2998 

1.2125 

545.2000 

0.4929 

683.0999 

0.5515 

695.8999 

1.0859 

743.0000 

0.9086 

1043.3999 

1.0897 

781.2998 

0.8292 

552.2000 

0.5066 

768.2998 

0.9325 

659.2000 

0.9111 

1282.2998 

1.2384 

767.7998 

0.6995 

609.2000 

0.4810 

681.7998 

1.0879 

725.3999 

0.8956 

925.7000 

0.8321 

523.8999 

0.5425 

704.2998 

0.7786 

562.2000 

0.4558 

694.5999 

0.8220 

934.2000 

1.1240 

1699.7000 

1.9828 

971.5000 

1.3375 

734.7998 

0.9847 

549.8999 

0.5234 

688.8999 

0.8782 

1369.2000 

2.0692 

849.7000 

1.1444 

558.0000 

0.4664 

551.2000 

0.5754 

743.0000 

0.7145 

1139.5000 

1.5717 

1042.5999 

1.3151 

941. 7000 

0.9080 

862.5000 

0.9340 

987.3999 

1.2192 

1085.2998 

1.1666 

728.3999 

0.8493 

1248.7998 

1.7349 

881.0000 

1.2301 

729.2998 

0.9287 

1062.0000 

0.8504 

1085.8999 

1.2003 

749.8999 

1.1255 

1136.5999 

1.1074 

579.0000 

0.7533 

998.2000 

1.1210 

625.0000 

0.8003 

642.5000 

0.8414 

713.0000 

0.9448 

625.5000 

0.6790 

967.5999 

1.2432 

1022.7000 

1.4037 

1006.7000 

1.3606 

596.7000 

0.6456 

1000.5999 

1.3256 

660.0999 

0.8086 

587.2998 

0.7184 

1051.5000 

1.2500 

658.8999 

0.8999 

711.3999 

0.8988 

677.5999 

1.0285 

1231.5000 

1.4286 

1004.3999 

1.335 

595.3999 

0.4473 

1255.7000 

1.8550 

800.7000 

1.1665 




TABLE  A.3(Contlnifd) 


ftinii 

LAMBDA 

1076.2000 

1.2229 

1277.7998 

1.8438 

1109.8999 

1.3950 

802.2998 

0.9201 

834.0999 

1.2543 

989.0999 

1.0521 

1392.0000 

1.9634 

1038.2996 

1.3092 

1250.7000 

1.4110 

829.3999 

0.9717 

728.2000 

0.8217 

807.3999 

1.1283 

1235.5999 

1.2762 

466.2996 

0.4801 

1047.2000 

1.2111 

1119.7998 

1.4488 

1199.2000 

1.6194 

725.8999 

0.8500 

652.7998 

0.7714 

868.5999 

1.3025 

723.8999 

0.7434 

734.7000 

0.7315 

681.2000 

0.8583 

1006.8999 

1.3533 

939.2000 

0.9735 
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